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Abstract. The dynamics of electric field generation and radial acceleration of ions 
by a laser pulse of rclativistic intensity propagating in an underdense plasma has been 
investigated using an one-dimensional electrostatic, pondcromotive model developed 
to interpret experimental measurements of electric fields [S. Kar et al, New J. Phys. 9, 
402 (2007)]. Ions are spatially focused at the edge of the charge-displacement channel, 
leading to hydrodynamical breaking, which in turns causes the heating of electrons and 
an "echo" effect in the electric field. The onset of complete electron depletion in the 
central region of the channel leads to a smooth transition to a "Coulomb explosion" 
regime and a saturation of ion acceleration. 
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1. Introduction 

Self-focusing and guiding of laser pulses is one of the most peculiar effects in 
nonlinear optics. In the case of superintense laser pulses (having irradiances beyond 
10 18 Wcm~ 2 /im 2 ) propagating in a plasma, self-focusing arises due to the combined 
effects of the intensity dependence of the index of refraction for relativistic electron 
velocities and of the expulsion of plasma from the propagation axis driven by the radial 
ponderomotive force [TJ El [3]; this latter effect leads to self-channeling, i.e. to the 
creation of a low-density channel that may act as a waveguide for the laser pulse. As 
the ponderomotive force acts on the plasma electrons, the channel is first drilled in 
the electron density only and becomes strongly charged [H [5] ; the electric field in the 
channel leads to ion acceleration [6j[71[8], which is of interest both as a diagnostic of the 
interaction and as a way to provide fast ions for applications, e.g. for the production of 
fusion neutrons [HI [101 [Hj . Besides these applications, the self-channeling process may 
be of relevance for laser-plasma acceleration of electron [12] and ion beams [T3l [H] , X- 
and 7-ray sources pJ2,[16], and fast ignition in laser-driven Inertial Confinement Fusion 

To interpret experimental results on self-channeling and related ion acceleration 
or neutron production, a simple modeling of the radial dynamics based on the 
ponderomotive force and electrostatic field only has often been used [TEJ El El [3, fTUl fT9]. 
Such a ponderomotive, electrostatic model (PEM) is attractive due to its simplicity 
and easy numerical implementation with respect to multi-dimensional, electromagnetic 
particle-in-cell (PIC) simulations which typically require massively parallel computing, 
while simulations based on the PEM in one spatial dimension (ID), taking only the 
radial dynamics into account, can be performed on a personal computer. Of course 
the simple PEM gives a strongly simplified description of the interaction where effects 
such as the nonlinear evolution of the pulse due to self-focusing or instabilities are not 
included. The use of the PEM can be justified only by its ability to reproduce, at least 
qualitatively, experimental observations for particular regimes. 

Recently, the proton imaging technique [20] allowed for the first time the detection, 
with spatial and temporal resolution, of the electrostatic, slowly-varying fields produced 
during and after the self-channeling process (previous experimental investigations 
[H EH E21 UHl (231 E] were based on optical diagnostics or indirect measurements not 
directly sensitive to electric fields). The proton imaging data were well reproduced 
numerically by simulating the proton probe deflection in the space- and time-dependent 
electric field distribution obtained from PIC simulations based on the ID PEM. 
Moreover, a similar dynamics of the electric field and the ion motion in radial direction 
was also observed in two-dimensional (2D) electromagnetic PIC simulations of the laser- 
plasma interaction [191 EI] • 

Motivated by the fair agreement with experimental results and more complex 
simulations, we have used the PEM to gain an insight of the radial dynamics during 
and after the self-channeling process. Although there may be much additional physics 
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at play that is not included in the simple PEM, a study based on the latter can be 
useful to discriminate purely ponderomotive and electrostatic effects from those due to 
other contributions, e.g. self-generated magnetic fields [22], anisotropy and polarization 
effects [25], hosing instabilities [261 [27], an d so on. 

It turns out that the dynamics contained in the ID PEM is already quite rich. We 
focus our attention on the regime where there is not a complete depletion of electrons in 
the channel and the electric field almost exactly balances the ponderomotive force (PF) 
locally. From the point of view of ion acceleration we call this the "ponderomotive" 
regime, since the force on ions is proportional to the PF, while we call "Coulomb 
explosion" the regime in which, due to complete electron depletion, the ions move under 
the action of their own space-charge field only. Our analysis shows that the transition 
between the two regimes occurs rather smoothly. 

A prominent effect observed in the simulation is the spatial "focusing" of the ions 
at the edge of the channel, where they form a very sharp spike of the ion density. The 
density spike splits up rapidly due to hydrodynamical breaking, and a short bunch of 
"fast" ions is generated. Density spiking and breaking occur even for pulse durations 
shorter than the time needed for ions to reach the breaking point. The onset of 
hydrodynamical breaking also causes a strong heating of electrons and the formation of 
an ambipolar sheath field around the breaking point. For pulse durations shorter than 
the time at which breaking occurs, there is a sort of "echo" effect in the electric field, 
which vanishes at the end of the laser pulse to re-appear later at the breaking location. 
Simple analytical descriptions of the spatial focusing mechanism and of the ambipolar 
field structure around the density spike are given. Finally, the smooth transition toward 
the "Coulomb explosion" regime is described both analytically and via simulations. 



2. The model 



We now give a detailed description of the ID, electrostatic, ponderomotive model which 
has been already used in Ref. [19] to simulate the radial electron and ion dynamics due 
to self-channeling of an intense laser pulse in an underdense plasma. Only the slowly- 
varying dynamics of the plasma electrons is taken into account, i.e. a temporal average 
over oscillations at the laser frequency is assumed. In other words, what we describe 
is the dynamics of electron "guiding centers", i.e. of quasi-particles moving under the 
action of the ponderomotive force (PF). The PF associated to a laser pulse described 
by the vector potential A(x, t) can be written, under suitable conditions (see e.g. |28j). 
as 

F p (x,i) = -m e c 2 V 7 (x,t), (1) 
7 (x,t) = [l + (a 2 (x,t)>], (2) 
a(x,t) = (e/m e c 2 )A(x,t), (3) 

where the angular brackets denote average over a period. We assume a non-evolving 
laser pulse (neglecting pulse diffraction, self-focusing and energy depletion) and complete 
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cylindrical symmetry around the propagation axis, taking only the radial dynamics into 
account. Under such simplifying assumptions, the laser pulse is completely defined by 
the cycle-averaged squared modulus of the vector potential in dimensionless units, which 
we write as 

a 2 (r,t) = (a 2 (x,t)> = a 2 {r)V{t) = a 2 ^ (4) 

where V(t) is the temporal envelope. The radial component of the PF can thus be 
written as 

F r = F r (r,t) = —m e c 2 d r ^, (5) 
7 = 7 (r,t)= [l + a 2 (r,t)] 1/2 . (6) 

Unless ao 1 the electron dynamics is relativistic. For ao •> 1, besides using the 
relativistic expression of the PF one has to account for the inertia due to the high- 
frequency quiver motion. This is included via an effective, position-dependent mass 
m = m e 7 of the quasi-particles [28]. We thus write for the radial momentum 

Pe, r - m e iv r . (7) 
Therefore, the equation of motion for the electrons is written as 

^ = F r - eE r (8) 
dt K 1 

where E r is the electrostatic field due to the space-charge displacement. The effect of 

the laser force on ions having mass = Am p 3> m e can be neglected, leaving the 

electrostatic force ZeE r as the only force on the ions. The ion equation of motion is 

thus written as 

^ = Z -^E r . (9) 

dt Am p v ' 

The electrostatic field is computed via Poisson's equation 

d r E r = A-Ke(Zn i -n e ). (10) 

Equations (El), (EJ) and (TIDl) are the basis of our particle simulations in ID cylindrical 
geometry. 

Thanks to the low dimensionality of our model, it is possible to use a very high 
resolution in the simulations. A typical run used 40000 spatial gridpoints, with a spatial 
resolution Ar = rf p /500 where d p = c/u> p is the plasma frequency, and up to 3 x 10 
particles for both electron and ion distributions. It turned out that such a high resolution 
is needed to resolve the very sharp spatial structures that are generated during the 
simulation, as it will be discussed below. The initial temperature of the plasma is taken 
to be zero and no significant numerical self-heating is observed during the simulations. 

In the simulation results shown below, the spatial coordinate r is normalized to the 
laser wavelength A, the time to the laser period T^, the density to the critical density 
n c , and the electric field to the "relativistic" threshold field E . The definition of these 
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parameters and their value in "practical" units for the typical choice A = 1 /im are as 
follows: 

T L = - = 3.34 x 1(T 15 s, (11) 

c 

2 

n c = = 1.11 x 10 21 cm" 3 , (12) 

Eo = 2tt^! = 3.21 x 10 10 V cm- 1 . (13) 
eA 

Momenta are given in units of m^c for ions and m e c for electrons, respectively. 

To compare with experiments, we note that the pulse intensity as a function of time 
and position is given by 

eA" 



I(r,t) = ^(A\r,t)) = J^-) a 2 (r,t). (14) 



JnA 2 N 1/2 



The relation 

a = 0.85 x ( - „ ) (15) 

V10 18 W cm- 2 /imV V ' 

thus gives the parameter CIq dbS Si function of the maximum intensity Iq, i.e. as the value 
of the intensity at the center of the laser spot and at the pulse peak. Some care should be 
taken when comparing to experiments where some average value, i.e. the pulse energy 
over the spot area and the pulse duration, is given instead; inserting such an average 
value into Eq. ffl5l) would represent too low a value for a proper comparison. 

The "pulse duration" and "spot radius" quoted in experimental papers most of the 
times refer to the FWHM of the temporal envelope and to half the FWHM of the radial 
profile of the intensity, respectively. In our model the laser pulse intensity varies as 

4 / nt \ ■ 4 ( L14t 



P( i) = sin*(-j=sin*(— y (16) 

for < t < t, while V(t) = for t > r. (A Gaussian pulse envelope has also been tested, 
but no significant differences in the simulation results was evidenced.) The parameter 
T1/2 = 1.14r/7r is thus the FWHM duration of the laser pulse intensity. For the Gaussian 
radial profile, the FWHM radius of the intensity profile is ri/2 = vhi 2r — 0.83r . 



3. Results 



The particle code based on the PEM was developed to analyze the experiment of 
Ref . [IH] , and the most relevant results emerged during such analysis. Therefore the first 
simulation we show, which gives us the basis for our discussion, has been performed for 
laser and plasma parameters in the range of those covered by the experiment of Ref.[19j. 
In the latter, the laser pulse wavelength was A = 1.053 /im, the intensity was in the 
range from 4 x 10 18 to 2 x 10 19 W cm -2 , the duration was n/2 — 1.2 ps ~ 330 Tl and the 
waist radius in vacuum was 7*1/2 — 5 A. As the effective values of the pulse amplitude 
and radius may change during the propagation into the plasma, the values of a and 
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Figure 1. (color online) Simulation results in the pondcromotivc regime. The profiles 
of electric field E r (blue, thick line) and ion density rij (red, dash-dotted line), and 
the phase space distributions of ions fi(r,p r ) and electrons f e {r,p r ) are shown at four 
instants (with times given in units of the laser period Tl) which are representative of 
subsequent stages of the dynamics: the initial ponderomotive acceleration (t = 38fTi), 
the vanishing of E r after the end of the laser pulse (t = 668Ti), the spiking of rij 
preceding the hydrodynamical breaking and the rebirth ( "echo" ) of E r correlated with 
strong electron heating (t = 8981l), and the "x" structure of the ion phase space 



after breaking (t 
T 1/a = 330T. 



1402Tl). Parameters are oq = 2.7, n e /n c = 10 , rp = 7.5A, 



ro were varied in the simulation until the best match in the reconstruction of proton 
images was found [19]. The plasma was created in a gas jet of Helium (charge state 
Z = 2, mass number A = 4) and the typical electron density n e was in the range from 
10 18 to 10 19 cm -3 , i.e from about 10~ 2 n c to 10~ 1 n c . 

3.1. Ion dynamics and hydrodynamical breaking 

FigIT] shows the spatial profiles of the electric field E r and of the ion density Hi, and 
the distribution functions in the (r,p r ) phase space of ions (fi) and electrons (/ e ) at 
four different times, which are representative of the subsequent stages of the dynamics. 
Figj2] shows the complete space-time evolution of E r and rii from the same simulation 
as contour plots. The parameters are a = 2.7, n e /n c = 0.01, r = 7.5A, Ti/ 2 = 330T L . 

The dynamics observed in the simulation of Fig{T] can be described as follows. 
During the laser pulse (t < r) , the PF E r pushes the electrons outwards, quickly creating 
a positively charged channel along the axis. This charge displacement creates a radial 
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Figure 2. (color online) Space-time contour plots of E r (r,t) (top frame) and rii(r,t) 
(middle frame) from the same simulation of FigfTJ The bottom frame shows enlarged 
view of the region marked by the dashed contour in the middle frame, showing the 
accumulation of ions on the axis at late times. The contour levels have been partly 
saturated to make regions of low field or density appear. 

ES field E r which holds the electrons back, as shown in the first frame of Fig{T] for 
t = 3817l (i.e. ~ 50 periods after the pulse peak). In this stage, we find the ES field 
to balance almost exactly the PF, i.e. eE r ~ E r (when plotting E r /e as well in FigJT] 
for t = 381, its profile cannot be distinguished from that of E r ). Thus, at any time 
the electrons are approximately in a state of mechanical equilibrium. From the electron 
phase space at t = 38 ITl we also observe that no significant electron heating occurs, f e 
being a narrow function along the During this stage the ions are accelerated 

by the electric force ZeE r = ZE r , and a depression in rii is thus produced around the 
axis. After the end of the pulse (t ~ 6687^ in Fig{T]), ion acceleration is over, the peak 
momentum of the ions is 0.015mjC, and E r ~ 0. This indicates that the electrons have 
rearranged their spatial distribution in order to restore the local charge neutrality. At 
the same time, we still observe a very weak heating of electrons, consistently with the 
keeping of the mechanical quasi-equilibrium condition up to this stage. 

However, the ions retain the velocity acquired during the acceleration stage. For 
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Figure 3. (color online) Detail of the profiles of electric field E r (blue, thick line) and 
ion density n, (red, dash-dotted line) around the "breaking" point at t = 898Tl, for the 
simulation of FigsQ] and [2l The peak value of rii exceeds forty times the background 
value no = 0.005n c . 



f > '"max, where r max ~ tq is the position of the PF maximum, the force on the ions 
decrease with r, and thus the ion final velocity does as well. As a consequence the ions 
starting at a position rj(0) > r max are ballistically focused towards a narrow region at 
the edge of the intensity profile. This spatial focusing effect is actually very tight, as 
in FigJT] we observe a large number of ions to reach approximately the same position 
(r = r& ~ 15.5 A) at the same time (t = tb ~ 8987^). Here the ions pile up producing 
a very sharp peak of rij. The peak value of rij is ~ 46no = 0.23n c at t = 8987^, as 
shown in the detail of the density and field profiles of Figj3] The density peak is well 
out of scale of the density axis in the plots of FigJTJ Simple modeling (see Appendix A ) 
provides an estimate of the position and the instant tb at which the density spike is 
formed clS cL function of the laser pulse parameters: 

/3\ 3/2 

n^i-j r ~1.8Ar , (17) 



h - -e^ «L ^ itfTj^. (18) 

2 V Z m e aoc V Z ao\ 

The latter expression is likely to be an underestimate for tb as the simple model assumes 

the pulse intensity to be constant during the time the ions take to reach the point r = r&, 

while in the simulation the pulse is shorter than tb and the value of the amplitude 

averaged over the pulse duration is lower than clq. For the run of Figfj] we obtain 

rb — 14A and tb — 569 Tl, in fair agreement with the simulation results. 

The piling up of ions at the point r = rb leads to hydrodynamical breaking in the 

ion density profile, as the fastest ions overturn the slowest ones at r = 7v The onset of 

hydrodynamical breaking is clearly evident in the contours of fi at the "breaking" time 
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t b ~ 898T £ (where the typical "vertical" shape of the contours of /j can be noticed). 
At later times (t = 1402T^) the ion phase space plot is reminiscent of the so-called 
"x-type" breaking that was first observed in the collapse of ion acoustic waves excited 
by Brillouin scattering [29j; in Ref.[2] the occurrence of this type of breaking at the walls 
of a self-focusing channel was also mentioned, but not discussed in detail. 

The ion density breaking leads to the generation of a short ion bunch (located 
near r = 22A at t = 1407Tj,), propagating in the outward direction, containing nearly 
monoenergetic ions having pi ~ 0.014mjC for the parameters of FigsHHS The motion 
of the bunch is almost purely ballistic, as it is evident in FigfSJ 

A similar feature was observed in the case of "longitudinal" ponderomotive 
acceleration described in Ref.[30], where the bunch formation is also shown to be related 
to spatial focusing of ions and hydrodynamical breaking of the ion density. In the case 
investigated in Ref.[30] ion acceleration is also of ponderomotive nature because the 
use of circularly polarized light prevents the generation of "fast" electrons and thus the 
onset of sheath ion acceleration. 

In addition to the ion bunch, we observe both a small fraction of the ions which 
is further accelerated at breaking (up to ~ 0.014mjc), and another fraction which is 
accelerated inwards, having negative pi up to ~ — 0.005mjC, and thus moves back 
towards the axis. At later times (t > 4000T^) these ions are found to form a local 
density maximum (i.e. a narrow plasma filament) around r = 0, as highlighted in FigfS] 
A local density maximum on the channel axis was also found in 3D electromagnetic 
simulations [22] but explained by magnetic pinching, which is absent in our model. 

The dynamics of the ions after breaking and related features, such as the "fast" ion 
bunch and the local density maximum on axis are related to the effect of the generation 
of a strong ambipolar field at the breaking point, which we now discuss. 

3.2. Electric "echo" effect and electron heating 

In the electric field plot at the breaking time we observe a strong ambipolar electric 
field appearing around the breaking point. The field is rather intense (its amplitude 
slightly exceeding that of the positive field due to charge depletion at earlier times) and 
highly transient; the complete "movie" of E r (r,t) in Figj2] shows that the field near the 
breaking point rises sharply from zero to its peak value over a few laser cycles time, 
and then decreases less rapidly to lower values (see the profile at t = 1402T^). The 
ambipolar structure slowly moves in the outward direction and is observable up to very 
long times. The "inversion" of the field, i.e. the appearance of a region where the electric 
fields points towards the axis, in the direction opposite to the initial stage of electron 
depletion, was evident in the proton imaging measurements reported in Ref.[19j. The 
electric "echo" is evidently correlated with the rapid and strong heating of electrons at 
breaking, which we observe in the f e plots at the breaking time. 

The generation of an electric field around the density peak as the ions approach the 
breaking point is interpreted as due to the inability of electrons to dynamically screen 




Figure 4. Electron energy distribution at r = 15. 8A and t = 8987l for the simulation 
of FigsQ] The dotted line is a two-temperature Maxwellian with Ti = 4 x lCP 3 m e c 2 
and T 2 = 2.5 x 10" 2 m e c 2 . 

the ion field when the density variation is too fast. Let d be the width of the ion peak at 
a certain time. The density changes in time due to the ballistic motion of ions moving 
towards the breaking point r = r&. The typical variation time of the ion density is 
Ti ~ d/vi, where u$ is the ion velocity determined by the acceleration stage. Electrons 
respond on a typical time of the order of the inverse plasma frequency, t p ~ 2ir/uj p , thus 
they are able to preserve plasma neutrality if d/vi > t p . If such condition is violated, 
the field around the density spike will be close to the unscreened field of the ions, i.e. a 
surface field of amplitude Ei ~ ±27rZeriid, leading to electron heating. The "neutrality 
condition" can be written in terms of the typical units used in the paper as 



(d/X) 2tt 




(19) 



which is useful to check its violation in the simulation results. However, it is noticeable 
that as a function of time d = d(t) may be not independent from the density n e = n e (t); 
due to mass conservation we expect d ~ l/n, and thus d ~ l/n e as far as n e ~ rij locally. 
Thus, as n e spikes the l.h.s. of Eq. (TT9ll decreases more rapidly than the r.h.s., boosting 
the violation of the inequality. This effect might account for the "robustness" of the 
electric "echo" effect in our simulations. The correlation between the formation of the 
density spike and electron heating might also be qualitatively understood on the basis of 
more general arguments, such as the fact that small oscillations of electrons around their 
point of (quasi-)equilibrium tend to become nonlinear across a sharp density gradient, 
i.e. around the spike, and lead to heating. 

The electron energy distribution f e near the breaking point at t — 898T^ (FigHj) 
for the simulation of of FigsH] may be roughly approximated by a two-temperatures 
Maxwellian with values Ti ~ 4 x 10~ 3 m e c 2 ~ 2 keV and T 2 ~ 2.5 x 10~ 2 m e c 2 ~ 12.8 keV. 
The value of Ti is fairly consistent with an electron acceleration in the unscreened ion 
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field leading to typical energies ~ eEid; in fact, estimating Ei ~ 0.05i?o and d ~ 0.1A 
from Figj3]we obtain a typical energy eE^d ~ 5 x 10~ 3 m e c 2 . The presence of the "hot" 
tail in f e may be a signature of stochastic processes able to accelerate a minor fraction 
of electrons to higher energies of the order of T 2 . 

The electron "temperature" generated during the highly transient heating stage 
accounts for the persistence of an ambipolar field around the density peak at later 
times. The almost Maxwellian shape of the electron energy distribution function allows 
to describe the late electric field as a "Debye sheath" field generated around a thin 
foil of cold ions and thermal electrons in Boltzmann equilibrium. In such a model 
the foil is "thin" if the spatial extension (FWHM) of the sheath field, £ s , is larger 
than the extension of the density "cusp" d, which is indeed the case in Figj3] where 
£ s ~ 0.3A > d ~ 0.1A. The model is analytically solvable in planar geometry (see 



Appendix B) for a delta-like foil density, which is appropriate in the limit d <C A D 



where A D is the Debye length. Using this analytical result for a rough estimate, we 
write the peak field E s and the sheath extension £ s as 

4A 2 

E s = 27renid, l s = 2—^, (20) 

d 

where is the ion density in the foil. Assuming parameters values from the 
ID simulations, nf ~ 0.23ra c , d ~ 0.1 A, we find in normalized units E s /E = 
7r(n/ /n c )(d/ A) ~ 0.07, fairly consistent with the simulation results. However, using 
T e = Ti ~ 4 x 10- 3 m e c 2 , we obtain \ 2 D /d = (T e /m e c 2 )(n c /n f ) (A 2 /47r 2 d) ~ 4.3 x 10~ 3 A 
and thus £ S /X — 0.034, which is smaller than the observed value and not consistent 
with the assumption £ s > d. This suggests that the sheath is mostly formed by the 
"hot" electrons having temperature T 2 ~ 67\ and lower density (roughly estimated to 
be n h — 0.1 7rif from the double-Maxwellian fit of the electron distribution), so that 
the effective values of \ 2 D and £ s may increase by one order of magnitude. Replacing rif 
by nh accounts for the screening by the colder electrons of the ion field acting on the 
"hot" ones. Due to these effects and additional ones (e.g. the dependence of the sheath 
profile on the cut-off energy for "truncated" Maxwellians, see Appendix B ) the simple 
thin sheath model cannot accurately describe the field around the density spike, but it 
is at least useful for a qualitative description. 

For clarity it is worth to point out that the breaking of the ion fluid and the 
formation of a strong ambipolar field around the breaking point occurs also for longer 
pulses, i.e. when the PF and the related electric field are not over at the breaking 
time. This has been observed by varying the pulse duration in our simulations. The 
electric field "echo" occurring when the pulse is already over is a remarkable signature 
of the strong electron heating that occurs following the ion fluid breaking in a regime of 
plasma neutrality. On a qualitative basis, a possible explanation of the electron heating 
is that near breaking the temporal variation of the plasma density becomes faster and 
a very strong density gradient is generated, so that electron oscillations around their 
equilibrium positions may become non-adiabatic, leading to heating. 

The generation of the strong ambipolar field structure, creating a very sharp field 
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gradient at the breaking point, has a feedback effect on the ion distribution, acting as a 
sort of "axe" that separates slower ions from faster ones. Slower ions are reflected by the 
negative part of the field and acquire negative radial velocity, thus driving the formation 
of the density maximum at r = at t > 1700 in Fig 121 Faster ions have enough energy 
to cross the electric field barrier, producing the escaping bunch, and may also receive 
additional energy from the positive part of the field; this accounts for the ions which are 
observed to get further momentum after breaking in Figfj] at t — 1402T^. These two 
ions components produce the upper right and lower left "arms" in the x-structure of the 
phase space at breaking. As the electric field decreases very rapidly after breaking, the 
other ions form a dense front moving in radial direction with a velocity much smaller 
than the fast ions in the bunch, as shown in FigfSJ Most of warm electrons remain 
confined around the ion front and their temperature is found to decrease with time. 



3.3. Transition to Coulomb explosion 

As stated above, under the action of the PF, the electrons are pushed away from the 
axis creating a back-holding electrostatic field which balances the PF almost exactly in 
the ponderomotive regime. However, the balance is possible only if the PF does not 
exceed the maximum possible value of the electrostatic force at some radius r, which 
occurs if all electrons have been removed up to such value of r, i.e. if complete electron 
depletion occurs, and ions have not moved significantly yet. Within our ID cylindrical 
model, an approximate threshold condition for complete electron depletion can be thus 
derived as follows. If all electrons are removed from a central region and the ion density 
is equal to its initial value, the electric field in the depletion region is given by 

Ed{r) = 27rZeriior = 2iren e or. (21) 

If F r exceeds the force due to the "depletion" field E d (r), this will first occur near r = 0, 
where F r is given approximately by 

F ~ a l (22) 

Hence, posing F r > eE d for r — > yields the condition 



2 

rl (l + al)^ >27re2n ^ (23) 
which we rewrite as 

1/2 

(24) 



a > 



k 2 / jU \ 1/2 



2 V 4 

2 



e 2 



m e c 



2n n e r = 2irr c n e r , (25) 



where r c = 2.82 x 10 13 cm. For n e = 10 19 cm 3 (10 18 cm 3 ) and ro = 7.5 /xm, we 
find k ~ 10 (1) and thus a complete electron depletion is expected to occur for ao > 10 



1.3]. 
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Figure 5. Threshold for the transition from the "ponderomotive" to the "Coulomb 
explosion" regime. The thick line gives amplitude threshold for complete electron 
density depletion in the central region [Eq. ([24|) ] as a function of the parameter n e o?"o 
with n e o in units of 10 18 cm -3 and ro in units of 1 /*m. Filled, empty and "gray" dots 
represent simulations where the electron depletion is absent, strong, or "marginal" , 
respectively. The labels a and b indicate two simulation for almost identical parameters 
but different pulse duration, so that in the longer pulse case (b) electron depletion is 
favored. 



The onset of complete electron charge depletion near the axis may occur even lower 
intensities than predicted by Eq. (l24"l) if the pulse is not too short. In fact, even if initially 
F r = eE r holds, as the ions move under the action of the force ZeE r = ZF r the density 
near the axis must decrease, so the maximum electrostatic field that can be generated 
near the axis also decreases in time, and the condition of complete electron depletion 
may be met. 

Figj5] shows the threshold amplitude as a function of n e r\. The data points in Figj5] 
represent the values of n e r\ for numerical simulations where electron depletion in the 
central region is either absent (empty, "white" diamonds) or evident (filled, "black" 
diamonds). With "gray" diamonds we indicate "near-threshold" cases where depletion 
occurs just very near to the r = axis and the maximum force on ions is still given by 
ZF r , i.e. the PF is larger than E r over most of the spatial range. The distribution of 
data points confirm that Eq. (1241) gives just a rough condition for the transition from the 
ponderomotive to the Coulomb explosion regime and that this transition also depends 
on the pulse duration. For example, the labels a and b in Figj5] indicate two simulations 
which have very similar values of ao and n e r\^ but for b the electron depletion is much 
stronger. This can be explained by the pulse duration for case b that is roughly two 
times the value for case a. 

In the Coulomb explosion regime, the ions in the region of electron depletion will 
be accelerated by their own space charge field. Figj6] shows results of a ID simulation 
where the laser pulse parameters are the same of FigH], but the density has been lowered 
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Figure 6. (color online) Simulation results in the regime of strong electron depiction 
and Coulomb explosion. All the fields and distributions are the same of Figfl] except 
that the normalized ponderomotive force F r /e is shown (blue, dotted line) in the 
leftmost upper plot. The vertical scales for the E r and pi^ r (ion momentum) are the 
same of Figfl] while to improve readability in the upper two plots for t = 902Ti and 
t = 1402Tl, E r has been multiplied by 10, and the scale for p e ^ r has been shrunken 
with respect to Fig[T] Parameters are ao = 2.7, n e /n c = 10~ 3 , ro = 7.5A, Tin — 330T. 



by a factor of 10 (thus, n e /n c = 10~ 3 ). For such parameters, from Eq. (1241) we expect 
to enter the Coulomb explosion regime. The plot of E r and of the scaled PF F r /e in 
Figj6] at t — 3807l and the corresponding contour plot of f e show that complete charge 
depletion occurs in the region r < r^ — 8 X. For r > r^, eE r ~ F r still holds. We notice 
that at the boundary of the depletion layer (r ~ r^), where the force balance breaks 
down, a few electrons are accelerated to relatively high velocity and escape towards the 
outer region. 

Due to electron depletion the resulting maximum force on the ions F max ~ ZeE r {r^) 
is less than the maximum of ZF r ; thus, the maximum momentum of ions at the end of 
the laser pulse (t ~ 7187^) is lower than in the case of Fig{T] for the same value of the 
laser intensity. Hence, for a given laser pulse, ion acceleration saturates for decreasing 
plasma densities as the condition ( }24"1) is crossed. The plots corresponding to the later 
times show that the spiking and breaking of the ion density, followed by ion bunch 
formation and ambipolar field generation, still occurs for these parameters; however, 
both the field amplitude and the electron energy are lower now. In simulations where 
the electron density is further lowered, electron depletion occurs over almost all the 
pulse profile, almost all ions are accelerated by the Coulomb explosion, while bunch 
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formation and ambipolar field generation tend to disappear. 
3.4- Saturation mechanisms for ion acceleration 

In the ponderomotive regime as defined above, the force on the ions is given by 
ZeE r = —Zm e c 2 d r [l + a 2 (r, t)] 1 ' 2 . If the temporal dependence of the pulse intensity 
is neglected, i.e. a 2 (r,t) ~ ct 2 (^), one may use the "ponderomotive potential", 
$ p (r) = Zm e c 2 [(l + a 2 (r)) 1 / 2 — 1] to estimate the final energy of ions as a function 
of their starting position. This assumption leads to a maximum ion energy equal to the 
peak value of $ p , 

11™ ^ $ P (r = 0) ~ Zm e c 2 [(1 + a 2 ) 1/2 - l] , (26) 

that is the energy acquired by ions initially located on the channel axis, i.e. at r = 0. In 
Ref.jH] this relation has been actually used to estimate the peak intensity on axis from 
the analysis of the ion spectrum. 

A clear limitation of this simple picture is the quite long time it would take for 
an ion starting at r = to go downhill the ponderomotive potential and to gain the 
maximum energy, as the PF is very small near r = 0. Thus, the laser pulse duration 
may be shorter than the acceleration time, limiting the ion energy. This is indeed the 
case for the simulation of FigJT] where the laser pulse is over before the fastest ions have 
moved out of the radial extension of the pulse. The maximum ion energy observed in 
Fig{T]is about four times lower than the estimate based on Eq. (126]) . In this case using 
the ion energy cut-off to evaluate a$ in Eq. (1261) would lead to an underestimate of the 
peak intensity. 

As discussed in section 13.3] the possible onset of complete electron depletion is 
another limiting factor for the ion energy as it gives an an upper limit to the accelerating 
force. The overall energy spectrum may also be modified by the onset of hydrodynamical 
breaking at the edge of the beam radius as discussed in sections 13.11 and 13.21 

To compare with preceding work we have simulated a case with almost the same 
parameters as those used in the calculations performed in Ref.jH] to interpret the 
experimental data. A He plasma is considered and the parameters are n e /n c = 0.01, 
7"i/ 2 ~ 1.2r = 6A, 71/2 = 300T, and a = 3.5 [this is lower than the value a = 5 quoted in 
Ref . [8] due to a different convention in the expression of the PF, i.e. F r = (1 + a 2 /2) 1 / 2 ]. 
Simulation results are shown in FigJT] The maximum ion energy, corresponding to the 
ions in the "fast" bunch, is ~ 2.2 MeV, very close to the energy spectrum cut-off in 
Fig. 2 of Ref. [8] , and significantly lower than the 2.7 MeV value corresponding to Eq. (]26]) . 
This is due to the complete electron depletion in the central region, occurring about 60 
cycles before the pulse peak and keeping the electrical force eE r below the maximum of 
the PF F r , as shown in FigfTJ The breaking of the ion density profile occurs at r ~ 11 A, 
very close to the analytical estimate given by Eq. ( 1A.4I) . and at t ~ 456Xl, roughly two 
times the value obtained from Eq. (lA.5l) . 
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Figure 7. (color online) Simulation results for parameters similar to those of Ref. [8] 
(ao = 3.5, n e /n c = 0.01, ro = 5A, t\ii = 300T). All the fields and distributions are 
the same of Figs 0] and [H 



4. Discussion and conclusions 

The one-dimensional electrostatic, ponderomotive model used in this paper to 
investigate the dynamics of self-channeling yields results whose agreement with 
experiments is remarkable, taking the simplicity of the model into account. A prominent 
example has been provided by proton imaging data which have been reproduced by 
simulating the particle deflection in the electric field computed by the present model 
|19j . Ion spectra have been also calculated for different laser and plasma parameters 
and agree with measurements within experimental uncertainties. 

Of course one should not forget that the model gives a very simplified description 
of the laser-plasma interaction, neglecting effects such as pulse diffraction and nonlinear 
evolution, and so on. The cases in which the simple model is successful in reproducing 
quantitatively or, at least, qualitatively some features observed in experiments or in 
more complex and self-consistent simulations, correspond to regimes in which the 
plasma dynamics during and after the self-channeling of the laser pulse is dominated 
by ponderomotive and electrostatic forces, other effects playing a secondary role. An 
example has been provided by the 2D electromagnetic simulations reported in Ref. [24] 
where the breaking of the channel walls, which has been characterized in detail with 
the ID model, causes the formation of two secondary laser beamlets propagating at an 
angle ~ A/r , consistently with a "leaking waveguide" picture. 

The analysis of the simulation results has evidenced details of the dynamics of 
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ion acceleration and electric field generation. In particular, hydrodynamical breaking 
has been shown to play an important role, causing electron heating, formation of 
an ambipolar field around the density cusp and, finally, affecting the ion spectrum. 
Unfolding this dynamics provided an insight on the formation of "x-type" structures in 
the ion phase space which had been previously observed in different contexts [29]. The 
related production of a dense, quasi-monoenergetic "bunch" of ions revealed similarities 
with the process of ponderomotive acceleration of ions in overdense plasmas [30J. 
A prominent consequence of hydrodynamical breaking, occurring for pulse durations 
shorter than the breaking time, is the "echo" observed in the electric field, i.e. its 
sudden rebirth after having disappeared at the end of the laser pulse. 

Throughout our work we found useful to distinguish the regime of "ponderomotive 
acceleration" from that of "Coulomb explosion" . In the first case, complete depletion of 
the electron density does not occur, the ponderomotive and the electric forces balance 
almost exactly and electrons are in a state close to mechanical equilibrium at all times 
before breaking; in such a regime the ponderomotive potential can be used to estimate 
the ion energy given the laser intensity (or vice versa), although the effects of finite 
pulse duration must be considered. In the second case, electron density depletion occur 
near the axis and the ions in this region are accelerated by their own space-charge field, 
leading to a saturation of the peak ion energy versus the laser intensity. An approximate 
analytical criterion for the transition between the two regimes have been given and tested 
by simulations, which also shows that such transition occurs smoothly. 

Acknowledgments 

This work has been supported by the Royal Society (UK) via a Joint Project, by the 
Ministry of University and Research (Italy) via a PRIN project and by CNR (Italy) via a 
RSTL project. A. M. acknowledges Queen's University, Belfast, UK for an International 
Fellowship. We are grateful to Alessandra Bigongiari, Peter Mulser, Francesco Pegoraro 
and Vladimir T. Tikhonchuk for useful comments and discussions. 

Appendix A. Analytical estimate of breaking point location 

A simple model can be used to account for the spatial focusing effect and estimate its 
characteristic parameters (i.e. the position and the instant at which the density spike is 
formed). For the sake of simplicity, let us neglect the temporal variation of F r . If F r was 
a linear function of r, the ion equation of motion would be of the harmonic oscillator 

type, 

m i d ^ = ZF r ~-k{r i -r b ), (A.l) 

implying that all the ions starting from an arbitrary radius rj(0) < r& would get to r = 
at the time tj, = (27r/fl)/4, where f2 = yjk/rrii. In our case F r is not a linear function 
of r, but a linear approximation of F r is quite accurate around the point r = r/ such 
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that d 2 F(r)\ r=rf = 0. We thus estimate the parameters r b and k in Eq. flA.ip from such 
a linear approximation. To further simplify the derivation, we take the non-relativistic 
approximation and write 

F r ~ Zm e c 2 d r (a 2 (r)/2), a(r) = a^ 12 ^. (A.2) 

The non-relativistic approximation turns out not to be very bad even if a — 1 because 
at r = rf > To the exponential factor is already small. By differentiating F r for two 
times, we easily obtain r/ = ^3/2ro. Expanding in Taylor's series around r = r/ we 
obtain 

k = F'(r f ) = Zm e c2 fo e -3/2 (A3) 
mi r 2 

Fir 13 /3\ 3 ^ 2 
r b = r f + -\H = -r f = - r ~ 1.84r . (A.4) 



k 2 1 \1 / 

Thus, r b depends on r only, and for the run of Figfj] we obtain r b — 14A. From the 
breaking time we obtain 



4 = - A /- = -e 3 /V4^. (A.5) 
2 V Ttj 2 V Z m e age 

This is likely to be an underestimate for the breaking time, since in our case the pulse 

duration is shorter than the ion acceleration time and it may be proper to replace the 

peak value of the amplitude Oq with some time-averaged value a < a . 

The predictions of the rough linear approximation of F r are thus in fair agreement 
with the simulation results. The important point to stress is that in the ponderomotive 
regime r& depends on ro only, while tf, depend only on the laser pulse parameters and 
on the Z/A ratio, but not on the plasma density. Our numerical simulations performed 
for different parameters in this regime show that the spatial focusing and piling-up of 
ions is a robust phenomenon which, once the spatial form of F r is fixed, tends to occur 
always at the same point. 

It is worth to notice that these estimates have been obtained for a Gaussian intensity 
profile. A different functional form would produce different results for r& and tj,, but 
their scaling with pulse width and amplitude should be the same. In general we expect 
any reasonable "bell-shaped" profile of the laser pulse to produce a sharp density 
increase near the edge of the beam, as this is the result of the decreasing with radius 
of the ponderomotive force in such region. The spiking of the density is sharp for 
ponderomotive force profiles such that d 2 F(r) = at some point. 



Appendix B. The sheath field around a thin foil 

In this section we compute analytically the sheath field around a plasma foil having a 
thickness much less than the Debye length, using a one-dimensional, cartesian geometry. 
The plasma foil is modeled as a delta-like distribution with the ion density given by 

rii(x) = nod5(x), (B.l) 
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where n Q d is the surface density of the foil. Ions are supposed to be fixed. The electrons 
are assumed to be in Boltzmann equilibrium 

n e = n exp ( -— ) = n exp ( ^ ] , (B.2) 



where $ is the electrostatic potential. It is convenient to express the coordinate, the 
potential and the electric field in dimensionless form, 

Z =^~, <P=7fT, £ = ^T' ( B - 3 ) 

A D J-e J-e^n 

where A 2 = T e / (47re 2 rio) is the Debye length. From Poisson's equation we thus obtain 

0" = e+. (B.4) 

We expect a potential that will be a even function of z (so that n e is even and E is 
odd), thus we can restrict the analysis to the z > region. 
Multiplying by 0' and integrating once we obtain 

i[0'(z)] 2 -i[0'(O+)] 2 = e ^)-e^). (B.5) 

Here 0'(O + ) = — Eq is the electric field at the surface. If the system is neutral, the field 
at z = +oo is zero, thus 

_ leg = -i[0'(O+)] 2 = e^ + °°) - e*°>. (B.6) 

Our first ansatz is to take e^ +OQ ^ = 0, from which we obtain [0'(O + )] 2 /2 — e^ = 0, and 
Eq. flB.5l) becomes 

l -[<P\z)f = e^\ (B.7) 

This latter equation can be integrated by the substitution = ln(/), i.e. f(z) = e^ z \ 
We obtain the potential 

= In ' 

and the electric field 



(z + 2/e y\ 



-21n(z/2 + l/e )), (B.8) 



e(z) = -0' = JJ^JJ^- (B.9) 
From Gauss's theorem we have E = 27ien d. Noting that e = we finally find 

<z) = -0' = . 2 ... , (B.10) 
\z + 4/d\ 

where d = d/X D is the foil thickness in units of the Debye length. This solution is shown 
in FigEH 

In dimensional units and for the whole sheath 

E{X) = |x|+4A 2 // lgn(x) = l + \x\d/4\l Slgn{x) - (B ' n) 
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Figure Bl. The field profile in the case d/X D = 0.1 for a cut-off energy u = oo (thick 
line), u = 5 (dashed), and u = 8 (dashed-dotted). The inner plot shows the profiles 
for small positve values of z including the u = 2 case (dotted). 

It is worth to notice that in this solution the peak field, Eq = E(x = O 1 * 1 ), does not 
depend on T, and that the spatial extension of the sheath (i.e. the distance at which 
the electric field falls by a factor 1/2), given by 
4A 2 

L = ^, (B.12) 

a 

is inversely proportional to d: the thinner the foil, the larger the sheath. Our approach 
is valid if L ^> d, i.e. A D d. 

In the above solution, the fact that the field extends up to infinity is related to the 
energy spectrum of electrons extending up to infinite energies. In some situation it may 
be more physically meaningful to assume that the electrons have a maximum energy 
t^max = uT e and "truncate" the maxwellian up to such cut-off energy. Pursuing this 
latter approach, we restart from Eq. (IB.5p . The maximum of the potential energy will 
be equal to U max = uT e , i.e. the maximum of <fi is — u, and <fi will be constant beyond 
that point. Thus, Eq.( 1B.5l) now reads 

-\el = -\W^ + )f = e- u -e^\ (B.13) 
and thus the equation for becomes 



0' = _^2(e^-e- u ). (B.14) 
Using the usual substitution / = we obtain after some algebra 

-u + In {I + tan 2 



-u/2 

-z — arctan(— ^=e l 



e 1 , /£o iU/2 N 



(B.15) 



. y/2 y/2 

The cut-off occurs at the point z r where <ft = —u, i.e. when the argument of the tan 2 
function equals zero; we thus find 

z r = V2e u/2 arctan(i|e n/2 ). (B.16) 
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The electric field is given by 



e[z 



) = v / 2e 



-u/2 



tan 



-u/2 



V2 



arctanf— %e u ^ 2 ) 

y/2 



(B.17) 



and e(z r ) = 0. The boundary conditions at z = remain the same, thus Eq = d/2X u . 

Fig JBll shows the field profile for different values of the cut-off energy u. The sheath 
extension z r is strongly dependent on u. 



References 



[1] Guo-Zheng Sun, Edward Ott, Y. C. Lee, and Parvez Guzdar. Self-focusing of short intense pulses 

in plasmas. Phys. Fluids, 30:526-532, 1987. 
[2] W. B. Mori, C. Joshi, J. M. Dawson, D. W. Forslund, and J. M. Kindel. Evolution of self-focusing 

of intense electromagnetic waves in plasma. Phys. Rev. Lett., 60:1298-1301, 1988. 
[3] E. Esarey, P. Sprangle, J. Krall, and A. Ting. Self-focusing and guiding of short laser pulses in 

ionizing gases and plasmas. IEEE J. Quant. Elec, 33:1879-1914, 1997. 
[4] A. B. Borisov, A. V. Borovskiy, V. V. Korobkin, A. M. Prokhorov, O. B. Shiryaev, X. M. Shi, 

T. S. Luk, A. McPhcrson, J. C. Solem, K. Boyer, and C. K. Rhodes. Observation of rclativistic 

and charge-displacement self-channeling of intense subpicosccond ultraviolet (248 nm) radiation 

in plasmas. Phys. Rev. Lett, 68:2309-2312, 1992. 
[5] A. B. Borisov, A. V. Borovskiy, O. B. Shiryaev, V. V. Korobkin, A. M. Prokhorov, J. C. Solem, 

T. S. Luk, K. Boyer, and C. K. Rhodes. Rclativistic and charge-displacement self-channeling of 

intense ultrashort laser pulses in plasmas. Phys. Rev. A, 45:5830-5845, 1992. 
[6] G. S. Sarkisov, V. Yu. Bychenkov, V. T. Tikhonchuk, A. Maksimchuk, S.-Y. Chen, R. Wagner, 

G. Mourou, and D. Umstadter. Observation of the plasma channel dynamics and Coulomb 

explosion in the interaction of a high-intensity laser pulse with a He gas jet. JETP Lett., 

66:830, 1997. 

[7] G. S. Sarkisov, V. Yu. Bychenkov, V. N. Novikov, V. T. Tikhonchuk, A. Maksimchuk, S.-Y. Chen, 
R. Wagner, G. Mourou, and D. Umstadter. Self-focusing, channel formation, and high-energy 
ion generation in interaction of an intense short laser pulse with a he jet. Phys. Rev. E, 59:7042- 
7054, 1999. 

[8] K. Krushelnick, E. L. Clark, Z. Najmudin, M. Salvati, M. I. K. Santala, M. Tatarakis, A. E. Dangor, 
V. Malka, D. Ncely, R. Allott, and C. Danson. Multi-mev ion production from high-intensity 
laser interactions with underdense plasmas. Phys. Rev. Lett., 83:737-740, 1999. 
[9] G. Prctzlcr, A. Saemann, A. Pukhov, D. Rudolph, T. Schatz, U. Schramm, P. Thirolf, D. Habs, 
K. Eidmann, G. D. Tsakiris, J. Mcyer-ter Vehn, and K. J. Witte. Neutron production by 200 
mj ultrashort laser pulses. Phys. Rev. E, 58:1165-1168, 1998. 

[10] V V Goloviznin and T J Schep. Production of direct fusion neutrons during ultra-intense laser- 
plasma interaction. J. Phys. D: Appl. Phys., 31:3243-3248, 1998. 

[11] S. Fritzlcr, Z. Najmudin, V. Malka, K. Krushelnick, C. Marie, B. Walton, M. S. Wei, R. J. Clarke, 
and A. E. Dangor. Ion heating and thermonuclear neutron production from high-intensity 
subpicosccond laser pulses interacting with underdense plasmas. Phys. Rev. Lett, 89:165004, 
2002. 

[12] V. Malka, J Faure, Y Glincc, and A F Lifschitz. Laser-plasma accelerators: a new tool for science 
and for society. Plasma Phys. Contr. Fusion, 47:B481-B490, 2005. 

[13] M. S. Wei, S. P. D. Mangles, Z. Najmudin, B. Walton, A. Gopal, M. Tatarakis, A. E. Dangor, 
E. L. Clark, R. G. Evans, S. Fritzlcr, R. J. Clarke, C. Hernandez- Gomez, D. Neely, W. Mori, 
M. Tzoufras, and K. Krushelnick. Ion acceleration by collisionless shocks in high-intensity- 
laser"underdense-plasma interaction. Phys. Rev. Lett., 93:155003, 2004. 

[14] L. Willingale, S. P. D. Mangles, P. M. Nilson, R. J. Clarke, A. E. Dangor, M. C. Kaluza, S. Karsch, 
K. L. Lancaster, W. B. Mori, Z. Najmudin, J. Schreiber, A. G. R. Thomas, M. S. Wei, and 



Electric field dynamics and ion acceleration 



22 



K. Krushelnick. Collimated multi-mev ion beams from high-intensity laser interactions with 

underdense plasma. Phys. Rev. Lett, 96:245002, 2006. 
N.H. Burnett and G.D. Enright. Population inversion in the recombination of optically-ionized 

plasmas. IEEE J. Quant. Elec., 26:1797-1808, 1990. 
Antoinc Rousse, Kim Ta Phuoc, Rahul Shah, Alexander Pukhov, Eric Lcfebvre, Victor Malka, 

Sergey Kiselev, Frederic Burgy, Jean-Philippe Rousseau, Donald Umstadtcr, and Daniclc 

Hulin. Production of a kev x-ray beam from synchrotron radiation in rclativistic laser-plasma 

interaction. Phys. Rev. Lett, 93:135005, 2004. 
Max Tabak, James Hammer, Michael E. Glinsky, William L. Krucr, Scott C. Wilks, John 

Woodworth, E. Michael Campbell, Michael D. Perry, and Rodney J. Mason. Ignition and 

high gain with ultrapowcrful lasers. Phys. Plasmas, 1:1626-1634, 1994. 
K. Krushelnick, A. Ting, C. I. Moore, H. R. Burris, E. Esarey, P. Sprangle, and M. Baine. Plasma 

channel formation and guiding during high intensity short pulse laser plasma experiments. Phys. 

Rev. Lett., 78:4047-4050, 1997. 
S. Kar, M. Borghesi, C. A. Cecchetti, L. Romagnani, F Ceccherini, T V Liseykina, A Macchi, 

R Jung, J Osterholz, O Willi, L A Gizzi, A Schiavi, M Galimberti, and R Heathcote. Dynamics 

of charge-displacement channeling in intense laser-plasma interactions. New J. Phys., 9:402, 

2007. 

M. Borghesi, A. Schiavi, D. H. Campbell, M. G. Haines, O. Willi, A. J. Mackinnon, P. Patel, 

M. Galimberti, and L. A. Gizzi. Proton imaging detection of transient electromagnetic fields in 

laser-plasma interactions. Rev. Sci. Instrum., 74:1688-1693, 2003. 
P. Monot, T. August e, P. Gibbon, F. Jakobcr, G. Mainfray, A. Dulieu, M. Louis- Jacquet, G. Malka, 

and J. L. Miquel. Experimental demonstration of relativistic self-channeling of a multiterawatt 

laser pulse in an underdense plasma. Phys. Rev. Lett., 74:2953-2956, 1995. 
M. Borghesi, A. J. MacKinnon, L. Barringer, R. Gaillard, L. A. Gizzi, C. Meyer, O. Willi, 

A. Pukhov, and J. Meyer-ter Vehn. Relativistic channeling of a picosecond laser pulse in a 

near-critical preformed plasma. Phys. Rev. Lett., 78:879-882, 1997. 
R. Fedosejevs, X. F. Wang, and G. D. Tsakiris. Onset of relativistic self-focusing in high density 

gas jet targets. Phys. Rev. E, 56:4615-4639, 1997. 
Andrea Macchi, Alessandra Bigongiari, Francesco Ceccherini, Fulvio Cornolti, Tatiana V Liseikina, 

Marco Borghesi, Satyabrata Kar, and Lorenzo Romagnani. Ion dynamics and coherent structure 

formation following laser pulse self-channeling. Plasma Phys. Contr. Fusion, 49:B71-B78, 2007. 
N. M. Naumova, S. V. Bulanov, K. Nishihara, T. Zh. Esirkepov, and F. Pegoraro. Polarization 

effects and anisotropy in three-dimensional relativistic self-focusing. Phys. Rev. E, 65:045402, 

2002. 

G. Shvets and J. S. Wurtele. Instabilities of short-pulse laser propagation through plasma channels. 

Phys. Rev. Lett., 73:3540-3543, 1994. 
N. M. Naumova, J. Koga, K. Nakajima, T. Tajima, T. Zh. Esirkepov, S. V. Bulanov, and 

F. Pegoraro. Polarization, hosing and long time evolution of relativistic laser pulses. Phys. 

Plasmas, 8:4149-4155, 2001. 
D. Bauer, P. Mulser, and W. H. Steeb. Relativistic ponderomotive force, uphill acceleration, and 

transition to chaos. Phys. Rev. Lett, 75:4622-4625, 1995. 
D. W. Forslund, J. M. Kindel, and E. L. Lindman. Plasma simulation studies of stimulated 

scattering processes in laser-irradiated plasmas. Phys. Fluids, 18:1017-1030, 1975. 
A. Macchi, F. Cattani, T. V. Liseykina, and F. Cornolti. Laser acceleration of ion bunches at the 

front surface of overdense plasmas. Phys. Rev. Lett., 94:165003, 2005. 



